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Summary. The Poincare-Dulac normal form of a given resonant system is 
in general non unique; given a specific normal form, one would like to further 
reduce it to a simplest normal form. In this note we give an algorithm, based 
on the Lie algebraic structure of the set of normal forms, to obtain this. The 
algorithm can be applied under some condition, non generic but often met in 
applications; when applicable, it only requires to solve linear equations, and is 
more powerful than the one proposed in previous work by the same author [Lett. 
Math. Phys. 42, 103-114; and Ann. I.H.P. 70, 461-514]. 

Introduction 

Poincare-Dulac normal forms (and their hamiltonian counterpart, Birkhoff-Gus- 
tavsson normal forms) are central to much of our understanding of local non- 
linear dynamics. 

We refer to Q for an exposition of the fundamental aspects - and to || 
for an in-depth exposition - of normal forms theory. Readable and reasonably 



sized introductions to the matter are contained e.g. in |19|, 24, |25|, |26|, [55J, |36[ ; 
see also [ pj6| . 

We want to discuss a point which has been investigated by several authors 
and which has an obvious interest both for theory and for concrete (including 
computational) applications, i.e. the possibility of further reducing a system 
which is already in normal form (NF in the following). 

It is indeed well known that the NF classification is in general redundant 
for resonant systems, and different NFs can be conjugated. This fact, and the 
possibility of obtaining more favourable NFs by taking advantage of the freedom 
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in the choice of the generators for the normalizing change of coordinates, were 
already remarked by Dulac in [p| . 

Among recent approaches to this problem, we recall ||, |5[ ^], ||, |l^, [l4|, pl| 
p8| , |29| , |33[ Many of them employ a similar set of ideas; among these, the 
use of the higher order analogues of the standard homological operator (already 
suggested in Q), and the Lie algebraic properties of the set of vector fields 
in normal form with respect to a given linear part (already suggested in Q). 
Each approach has also, of course, its own features, and leads to similar but not 
identical results. 

In |^0|, [H] I proposed an approach based on iteration of the standard 
Poincare normalization procedure (the reduced normal forms thus obtained 
were therefore called Poincare renormalized forms, or PRF for short); this is 
completely algorithmic and requires only to solve linear equations at each step. 

The PRF algorithm uses only the grading by homogeneity degree of the 
algebra Qa of polynomial vector fields in normal form with respect to a given 
linear part A; it is thus quite general, but fails to take advantage of the Lie 
algebraic structure of Qa- Not surprisingly, taking this Lie algebra structure 
into account can lead to a considerable simplification of the computations and 
of the resulting simplified normal form. 

The main purpose of this note is indeed to conjugate the PRF approach with 
Lie algebraic considerations, so to obtain a procedure which takes advantage of 
the Lie algebra structure of Q and keeps the computational simplicity of the 
PRF approach. Rather than aiming at the greater generality, we will focus on 
a structure which is non-generic, but relatively common in applications. 

The discussion given in this note should also clarify some points related to 
standard PRF approach (as we discuss in the lines below); it should also make 
clear how much the PRF approach is related to Broer's ideas and to previous 
work along the same lines. Actually, the set of ideas employed here is essentially 
the same (higher order homological operators, Lie algebras nitrations) as in 
those works; they are blended here in a form which is suitable for concrete 
computational implementation. 

It should be stressed that concrete computations in [£20|, |2l| were performed 
by using this the Lie algebraic structure, although this point was not sufficiently 
stressed there; this could have caused some confusion, which the present note 
should hopefully dissipate. To make things worse and add occasions for confu- 
sion, in [^(J |2l[] the term PRF was on the one hand precisely defined, but on 
the other end also used to indicate generically reduced normal forms obtained 
by use of a sequence of Poincare transformations, such as those based on the 
Lie algebra structure. 

This can cause confusion due to the following fact: if we consider the se- 
quence of Lie-Poincare transformations prescribed by the PRF procedure, we 
obtain some reduced normal form F\ (which is of course a PRF). If the se- 
quence of Lie-Poincare transformations is not taken in this order but according 
to a different scheme, e.g. to take advantage of the Lie algebra structure on Qa 
(as in considering the above mentioned example in p^, pl|), we obtain another 
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reduced normal form F2 . In general F2 ^ but moreover this F2 can fail to be 
a PRF according to the proper definition; this happened for the main example 
in |20| |2lJ and will also be the case for the example B considered in detail below. 

The paper is organized as follows: in section 1 we recall the main facts 
about Poincare-Dulac (standard) normal forms and fix notation; in section 2 
we discuss two further reduction scheme, i.e. the "generic" PRF procedure 
and a modification of it which makes use of the Lie algebraic structure of Qa 
and applies when this structure has favourable properties (this will be the LRF 
approach); in section 3 we briefly consider a three-dimensional example, with 
linear part corresponding to a center-focus, showing the simplicity of compu- 
tations required by the LRF approach; in section 4 we consider a case where 
the structure of Qa is not the optimal one for use of the LRF approach, but 
still allows for its use. Finally, in section 5 we analize in full detail a simple 
two-dimensional example, giving a completely explicit description of the renor- 
malizing transformations and renormalized forms up to order six, and compare 
the results obtained with the PRF and with the LRF approach. 

We also provide two appendices: in the first one we discuss Bruno's treat- 
ment of PRF and remark that his definition is not equivalent to the one given in 
[ po[ |2f| ] , so that the example he considers in |L2], [l3| - and which falls in the case 
considered in section 5 - does not apply to PRFs; in appendix B we briefly recall 
the main ideas put forward by Broer and Baider, and describe the Broer-Baidcr 
reduction procedure in the language employed in the present paper. 
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1 Standard normal forms 

Let us first collect some basic formulas about (standard) normal forms. We 
work in R™ with basis {ei, ...,e„} and coordinates (a; 1 , ...,x n ); and consider a 
vector field X in R™ having a zero in the origin. This is written in coordinates 
as X — p(x)di (here and below, di — d/dx l ); we expand the vector function 
f(x) in homogeneous terms as f(x) = Y^T=ofk( x ), witn fk(ax) = a k+1 f k (x). 
When considering coordinate expressions, we will write Vk for the set of vector 
fields with components homogenoeus of order (fc + 1) in the x. It is clear that 
under commutator we have [14, V m ] C T4 +m . 

We also consider the linearization Ax of / in the origin, given by A 1 ^ :— 
{df 1 /dx ) (0) . We write in general f(x) = Ax + F(x), where F(x) collects 
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nonlinear terms only. 

We denote by a := {Ai, A n } the spectrum of A. If the Ai satisfy some 
relation 

n 

:= {pt ■ A) = A Q (1) 

i=l 

where a = 1, ...,n, /ij are non-negative integers, and :— Y^i^i — we say 
that A is resonant. 

A vector of the form v = (x^ 1 ...x^ n )e a , with ^ and a as in (1), is called 
a resonant monomial vector, the linear span of resonant monomial vector is the 
linear space of resonant vectors. When some ambiguity could arise, we specify 
these are resonant with A. 

The Poincare-Dulac theorem affirms that it is possible to find a sequence of 
near-identity changes of coordinates such that in the new coordinates (denoted 
again by x) the vector field X is expressed as X — g l (x)di, with g(x) — Ax + 
G(x), and with the nonlinear term G(x) being the sum of resonant terms only. 
The change of coordinates defined by the whole sequence is in general only 
formal, but it is analytic (in some open neighbourhood of the origin) if zero lies 
out of the convex hull of a in the complex plane (Poincare condition) . 

It is actually also possible to express the changes of coordinates as the time- 
one action of vector fields H k = h\di (this has several theoretical and practical 
advantages j8, 17, 30, 37|); in this case we speak of Lie-Poincare transformations. 



When a nonlinear vector field X = [A^x + G l {x)}di in given coordinates 
satisfy the condition that the nonlinear part G l (x) is a resonant vector - with 
respect to the linear part A - we say that X is resonant (in the x coordinates). 

Notice that the resonance condition involves only the eigenvalues of A; if 
we decompose A into a semisimple and a nilpotent part (mutually commuting), 
A = A s + A n , then only A s enters in the definition of resonant vectors. 

If A is not semisimple, it is actually possible to refine ulteriorly the form 
to which G(x) can be taken. Denoting X& = A^x^di the vector field associated 
to the linear part of X (notice this is invariant under near-identity changes of 
coordinates), we denote by Cq the operator of commutation with Xa] that is, 
C (Y) = [X A ,Y]. Obviously, C Q : V k -> V k , 

We can then force Gk = G D Vk to be in a complementary space to the 
range of Co in each Vk- If we introduce a scalar product^] in Vk (and hence 
in V = Vo U Vi U ...), then we can force G & Ker(£g ). When the condition 
Xq e Ker(£j) is satisfied, we say that X is in Poincare-Dulac normal form (in 
the given coordinates). 

Notice that in general [AT^X^+J / 0, and thus the adjoint of the linear 
part of a vector field in normal form does not commute with the full vector field, 
but only with the nonlinear part. On the other side, if LY4, AT^+J = (which 
is equivalent to [A, A + ] — 0, i.e. to A being a normal matrix), then X^+ is a 
symmetry of X . 

1 Choosing the Bargmann scalar product Jl9| will guarantee that Cq(Y) = [X A + , X], where 
A+ is the adjoint of A. Another popular choice (standard scalar product) is described in [El. 
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Here we will assume that A is semisimple A = A s and has been taken to 
Jordan normal /orm^|, so that 

A = diag(Ai,...,A„) . (2) 

Thus, in the present discussion, vector fields in normal form and resonant vector 
fields will be the same. We also have that £ = -Co (provided we have chosen a 
suitable, e.g. the Bargmann, scalar product). 

Resonant vector fields are thus identified by the condition 

[X ,X] = (3) 

where Xq is the linear part of X. In order to emphasize that we are actu- 
ally dealing with coordinate expressions and not with geometrical objects, it is 
maybe better to deal with components of X in the x coordinates. It is easy to 
see that if Xf = f l (x)di and X g = g l (x)di, then [Xf, X g ] = Xh = h z (x)di, with 

h(x) = {f(x),g(x)} := (P{x) ■ d^x) ~ (g j (x) ■ a, )/*(*) . (4) 

The bracket {., .} is just expressing the commutator in terms of the components 
of vector fields in given coordinates. Relation (3) reads then 

{[(Df)(0)]x,f(x)} = 0. (3') 

The point is that once f(x) is in normal form, we can still consider near- 
identity Lie-Poincare changes of coordinates taking it into a different normal 
form f(x); in order to be guaranteed that f(x) is still in NF, i.e. (3') is satisfied 
by /, the generator h must be chosen to be itself resonant with A = (Df)(0). 
Indeed, the set Ga of vector fields commuting with Xa is obviously a Lie algebra. 

We recall that if h(x) — x + YlkLi h>k{x), the changes of coordinates given 
by the time-one flow of Xh maps X into X = e H Xe~ H . This can be explicitely 
computed by the classical Baker-Campbell-Haussdorff formula; with Q(X) := 
[H, X] , we have 

oo 

X = £(1/ S !)9 S (X) (5) 
which in terms of homogeneous components reads 

[m/k] 

f™=T,- n S (fm-sk) (6) 
s=0 

where [m/k] denotes the integer part of (m/k) and TL{f) = {h,f}. 

It is convenient, for further discussion, to define the higher order homo- 
logical operators (already considered by Takens Q) as Cu '■= {/fe, •}• We also 
denote, for later discussion, by G = C(A) the centralizer of A in the Lie algebra 
of n-dimensional matrices (with Lie operation the matrix commutator) ; a basis 
for this will be given by matrices {K\, Kd\ (if A ^ 0, we can choose K\ = A); 
notice that d < n. 

2 The theory of normal forms is well developed with out the se assumptions, but reduction 
encounters a number of substantial obstacles; see e.g. |32l p6| for normal forms in the case 
where A n ^ and/or A is not in Jordan normal form. 
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2 Two non-equivalent further reduction 
schemes: PRF and LRF 



Having denned normal forms and obtained the formulas for the change of vari- 
ables, we can wonder if one can change variables thus transforming a normal 
form into another one, more convenient in some respect (in particular, having 
a smaller number of nonlinear terms, or a smaller number of low order ones). 
The answer is obviously yes, as already noted by Dulac | fl8|| , in the form and 
within the limits implied by (5), (6) above. 

In recent works - as mentioned in the Introduction - I have proposed a 
general procedure (and algorithm, requiring to solve only linear equations) for 
such a reduction |21j| . Since this approach represents a direct extension of 
Poincare procedure, the normal form (in general not unique) obtained in this 
way has been called "Poincare renormalized form" (PRF), and we refer to it as 
the PRF approach. 

On the other hand, in some cases it is also possible to proceed in a different 
way: one can use the Lie algebraic structure of the set of resonant vectors. We 
will call the reduced normal forms obtained in this way, "Lie renormalized form" 
(LRF), and we refer to the procedure as the LRF approach. 

It should be mentioned that the relevance of Lie algebraic structures in 
normal forms theory was stressed (in his thesis) by Broer, who gave a very 
general and powerful reduction procedure (see appendix B). 

Our procedure will be less general and powerful, and apply only in favou- 
rable cases (see below); but it is simpler and - when applicable - its practical 
implementation in completely explicit computations is elementary, as we also 
show by example. 



2.1 The PRF approach 

I will now recall the basic aspects of the PRF approach, developed in [^0[ pl[ ; 
the reader is referred to these works and to [jl6| for further detail. 

Let W be the vector field (say already in standard NF, to avoid trivial 
steps) under consideration. Let us write it as W = X + J2T=i Wk, with X 
linear and Wk homogeneous of degree (k + 1). For the sake of simplicity^ we 
will assume Xq ^ 0. 

Take the first nonzero Wk, say W Po ; operating with transformations gen- 
erated by e VVfc ^ :— Kcr(£ ) n (successively for k = 1,2,...) we can 
eliminate all terms in [W po , W^], i.e. in the range of M po , defined as the 
restriction of the operator C po to Ker(£ )- 

3 Notice that the PRF can as well deal with cases where the linear part vanishes: now the 
standard homological operator is only the first in a chain of operators, and we can use the 
other ones for reduction. In practice, computations will be very difficult for A = 0, unless 
some other constraint (e.g. symmetry properties) reduces the set of allowed nonlinear vector 
fields; see for PRF analysis of symmetric systems with A = 0. 
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Let now W Pl be the first nonzero term among those with k > po in the 
normal form obtained after the above transformations. Using transformations 
generated by h k ^ S Wjp ■— Ker(£ ) n Ker(£ Po ) n Vk (successively for k — 
1, 2, ...), we can eliminate all terms in [W Pl , VU^ ], i.e. in the range of M Pl , the 
restriction of C Pl to Ker(£ ) H Ker(£ Po ). 

The process can obviously be continued indefinitely, until either all nonlin- 
ear terms of degree higher than p q are killed, or Ker(£ )nKer(£ Po )n...nKer(£ Pii ) 
is empty. 

Notice that in this way due to the restriction to kernels of lower order 
homological operators (i.e. due to the use of M. p rather than C p ), at each step 
we are not affecting the terms which have already been simplified; in facts at 
each step we stabilize new terms. 

The generators ni are chosen as solutions to higher homological equations: 
if /P^ is the term of order (A;+l) after the first j rounds of further normalization, 
and 7r["^ the projection from Vk to Wu \ this is given by 



[m/k] 

2J ~ T~L S (fm-sk) 

3=0 S - 



(7) 



0. 



where H — {h^ \ .}; that is, see (5), by requiring that (wjf^ 

(1) 

Notice that the h\ considered above are in general not uniquely defined: 
they are unique up to an element in Ker(A4 Pj ), see above. Thus, the PRF is in 
general not unique. 

The reduced normal form obtained according to this procedure will be 
called a Poincare renormalized form (PRF) for the normal form W = W^°\ A 
precise description of the spaces to which the terms non-eliminable in this way 
belong - i.e. an abstract general description of PRFs - and a more detailed 
discussion are given in (M E0| El] . 



2.2 The LRF approach 

Let us now consider a different further normalization scheme. 

Consider the set of vector fields in R n which are in normal form with respect 
to the given linear part A, i.e. the set of Y G V such that = 0. It is 

obvious that these form a Lie algebra (the Lie operation being the standard 
commutator of vector fields); we denote this algebra by Q. 

Let us recall a general characterization of vector fields in normal form rele- 
vant in this context [36|. Consider the linear vector field X^; we say 
that the differentiablc function ip : R™ — » R is an invariant for if X^((p) = 0. 

Denote by I* (A) the set of invariants for X& which are meromorphic (that 
is, can be expressed as a quotient of algebraic functions) in the x coordinates; 
denote by T(A) C 1* (A) the set of algebraic invariants for Xa, and by Tk(A) C 
1(A) the set of algebraic invariants for Xa which are functions homogeneous of 
degree k + 1 in the x variables. 
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Let G = C(A) be the centralizcr of A in the algebra of n dimensional 
matrices; let its Lie algebra be spanned by matrices {K\, ...,Kd} (we can always 
assume K\ = I, and that K a = A for some a, provided A =^ 0; notice that 
d < n). We denote by X^ the vector fields corresponding to these, i.e. given 
in the x coordinates by X 1 -^ — (K a x) l di. 

Then the most general vector field W in Q can be written as 

d 

w = M*) xia) (8) 

a=l 

where fi a (%) G I* (A). In other words, Q is contained in a finitely generated 
module over 1* (A) . 

Notice that the vector field W must be algebraic in the x, and X^ are 
linear in x, so that functions fi a (x) 6 I* (A) having poles of degree d > 2 in 
x = cannot appear in (10). That is, only algebraic functions and functions 
with simple poles in the origin can appear in the actual normal form unfolding: 
Q is not the full G-generated module over I* (A) . 

Example. Let us briefly mention an example where indeed meromorphic func- 
tions of the invariants enter in the normal form unfolding. Consider systems in 
R 3 with coordinates (x, y, z); let the linear part be given by the diagonal matrix 
A = diag(— 1, 1, 2), so that Xa = —xd x + yd y + 2zd z . This has two basic invari- 
ants, given by \&i = xy and \&2 = x 2 z. We take as X^ the vectors X^ = xd x , 
XW = ydy, and X^ = zd z . It is immediate to check that (*? /* 2 )^ (3) = y 2 d z 
and (^/2/^i)X^ = xzdy are polynomial and resonant with Xa- 

In several cases it happens that Q has a more convenient structure, i.e. the 
Ha in (10) can actually be taken to be in 1(A), and not just in I*(A). In this 
case we say that all the vector fields in Q are quasi-linear, or that we have a 
quasi-linear normal form. In particular, this is the case when A admits only 
one basic invariant (see the examples below). 

If the normal form is quasilinear, we have Q n Vk+i = Ik{A) <8> G, and 
the analysis of the structure of Q results to be particularly simple, as we now 
discuss. 

Call X* the algebra spanned by vectors which are written as X — s(x)X^ 
with s E I* (A); call X a the algebra spanned by vectors as above with s G 1(A) 
(this is the module over T(A) generated by X^). 

As seen before, in general we have X x © ... © X d C Q c X* © ... © X%, and 
in the (favourable) quasi- linear case we actually have Q = X\ © ... © X^. 

Consider now the commutation relations between elements of the subalge- 
bras X* and X^; it is immediate to check that 

[^)X^,a f ^)X^] = (^Wido-p/d^X^i^XW 

(9) 
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Notice that when X^ = X A , by definition JfW(^) = 0, and [X (a \ iW] = 0; 
thus the corresponding subalgebra Af^ is always an abelian ideal in Q. 

Remark. Note also that, as obvious from the formula (9) above, the union of 
subalgebras X ai U ... U X as is a subalgebra in Q if and only if {X^ ai ^ , ...,X^ a ">} 
span a subalgebra in G. 

It can happen that we are able to determine a sequence of subalgebras 
J~ p C Q, each of them being the union of X a subalgebras, such that Tq — Q and 

\Q,T P \ = T P+ i ; (10) 

if this terminates in zero we say that Q has a quasi-nilpotent structure. Notice 
that the factor algebras T p := J- p /J- p+ i are in general not abelian. 

By the above remark, Q can have a quasi-nilpotent structure only if G is 
nilpotent. The chain of subalgebras T v C Q can then be read off the descending 
central series G p of G; recall that the factor algebras 7 P = G p /G p +\ for this are 
abelian. The subalgebras T p introduced above are therefore moduli over 1(A) 
generated by abelian subalgebras j p of G. 

Assume now Q is quasi-nilpotent. In this case we can first work with gener- 
ators in ri and simplify terms in I\ (e.g. by following the PRF algorithm within 
the set this allows to work with more familiar projection and homological 
equations than if setting the problem in a completely Lie algebraic framework) , 
then consider generators in and simplify the corresponding terms being guar- 
anteed that Ti terms are not changed, and so on. 

Notice that in this case we are - roughly speaking - just using the nilpotent 
structure of (the finite dimensional group) G, rather than the one of (the infinite 
dimensional algebra) Q. 

Needless to say, this approach is particularly convenient when the T p are 
generated by a single element of G. 

The situation depicted above is met in applications: e.g., it applies to any 
nontrivial two-dimensional case and several three-dimensional ones p^| . More 
generally, it always applies when there is only one basic invariant. 

In the following sections we will consider some simple examples where the 
LRF is easily computed, and it turns out to be definitely simpler than the PRF. 



3 Example I 

Let us consider a system in R 3 (we use coordinates x, y, z) with linear part given 

by 




It is easy to see that this has only one basic invariant W := (a; 2 + y 2 ). As any 
meromorphic function of \E' is either algebraic or has a pole of degree d > 2 in 
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the origin, wc deduce that the most general vector field in normal form with 
respect to this linear part is 

oo 

W = X A + Y, a kX k + bkY k + c k Z k (11) 

k=l 

where (with k > 0) 

f X k := (xd x + yd v ) 
{ Y k := * k {-yd x +xd v ) 
{ Z k := ( zdz ) 

This form can also be easily deduced by explicit computation applying the 
definition of resonant vector field. Obviously Xa = Y a — Zq. We denote by 
H,v,cf the first k > 1 such that a k ,b k ,c k are nonzero. 
The X k , Y k , Zk satisfy the commutation relations 

[X k ,X m ] = 2(m - k)X k+m , [Y k ,Y m ] = , [Z k , Z m ] = 
[X k ,Y m ] = 2mY k+m , [X k ,Z m ] = 2mZ k+m , [Y k , Z m ] = 

Denoting by X , y, Z the algebras spanned by the X k , the Y k and the Z k , 
we have that Q = X © y © Z, and that y © Z is an abelian ideal in Q. We can 
thus apply the LRF procedure discussed above. 

We first operate on X, with generators also in X (thus H k = a k X k ): in 
this way we can eliminate all terms except the X^ and the X 2f _ t ones, as implied 
by (12). In doing this we modify terms in y © Z. 

Having performed this first step, we pass to consider the y and Z terms, 
operating with generators in y © Z (thus H k = j3 k Y k + "f k Z k ). It is clear from 
(12) that we can eliminate all terms with k > /U, but no lowest order ones. Thus 
we end up with a LRF given by 

_ H fJ. 

W = X A + 0,^X^ + 0,2^X2^ + ^b k Y k + y^c k Z k ; 

k— v k—cr 

the hat on constants mean that these are not the same as in the initial form 
(11), and obviously a sum with lower limit greater than the higher limit should 
just be meant as zero. 

As shown by this example, the computations required for the determination 
of the general LRF are actually very simple. 

4 Example II 

We want now to consider an example where the algebra Qa of resonant vector 
fields has not the optimal structure for LRF reduction, i.e. the T p are not 
generated by a single element of G (see section 2). 

We consider a system in R 4 with linear part corresponding to two oscillators 
with nonzero and non-resonant frequencies a and (3, \a/{3\ Q, i.e. in block 
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notation 

aJ \ T /0 -1 



(3 J J ' 3 \\ 

this A has eigenvalues A = ±ia, ±i@. We use coordinates (x, y, z, w), so that 

= a(-yd x +xd y ) + (3 (-wd z + zd w ) . 

It is immediate to see that the system is simply resonant and admits two inde- 
pendent basic invariants, tpi — x 2 + y 2 and tp 2 — z 2 + w 2 . On the other hand, 
the linear space Qi of linear vector fields commuting with Xa is spanned by 
X^ — (M Q x)V, with a = 1, ...,4 and we can choose the matrices M a e.g. as 
(in block notation, with / the two-dimensional identity matrix and J as above) 

»)- M -(o S)- M -(o I) ■"--(I ' 

It is also immediate to check that X^ a \i})i) — except for X^tyi) = 2tpi and 
X( 2 Hi; 2 ) = 2Tp 2 . 

We will write X^ m :— ipitp™ X^ ; these span 5, for fc, m non-negative 
integers. We also write, with the notation introduced above, Q = X\ ® ... ® X4; 
note that X% © A4 is an abelian subalgebra in G. 

The nontrivial commutation relations are given by 



vl x 1 

km' pq 

X 1 X 2 

km' pq 
yl va 
^km' ^pq 



- 2(p k)Xl +p m+q , [X| m ,Xp ? ] — 2(q m) X 2 +p m+q 
= 2p X k+pm+q — 2mX k+pm+q 

2pX£ +P , m+ q . [Xl n ,X°]=2qXZ +p>m+g (a = 3,4) 



It is easy to see that there is no way to separate the algebra as the sum of 
one-dimensional moduli over 1(A). Thus, the best we can do in the spirit of 
the LRF approach is to act first in the X\ © X 2 subalgebra. Once this has been 
reduced, say with lowest nonlinear term in Q^, we can eliminate all terms in 
(X 3 © X4) n G m for all m > fj,, as follows immediately from the commutation 
relations. 

Note that the choice of terms to be eliminated in X\ © X 2 is to a large 
extent arbitrary, and correspondingly the LRF is not unique. 



5 Example III 

We will now consider the LRF procedure for a simple system in full detail, i.e. 
aim at producing completely explicit formulas for the renormalizcd form and for 
the renormalizing transformation, thanks to the computational simplicity of the 
LRF procedure. We will also compare explicitely the LRF and PRF reductions. 

Let us consider a two-dimensional system (in R 2 with coordinates x, y) 
with linear part given by 



i.e. corresponding to the vector field yd y . We note immediately that here A is 
diagonal, so that A = A + . We have only one basic invariant tp(x,y) = x. 
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5.1 Normal forms 

It is easy to see that the kernel of Cq is spanned by the arrays of vector fields 
(with k > 0) 



(with this notation the linear part considered here is given by Yq). These vector 
fields satisfy the commutation relations 



[X k , X m ] = {m-k) X k+m , [X k , Y m ] = m Y m , [Y k , Y m ] = . (13) 



We denote by X the algebra spanned by the X' k s, by y the algebra spanned by 
the Yfc's; obviously Q = X ®y . Note that y is an abelian ideal in G. 

The (standard) normal form corresponding to the linear part considered in 
this section will thus be given by a vector field 



depending on the two infinite sequences of real constants a kl b k . 

We will denote by fi (respectively, by v) the first k > 1 such that a k ^ 
(respectively, such that b k ^ 0). It results that for the sake of our present 
discussion, it is of interest to consider the case v < fj, (for a full discussion of 
this system, whatever [i and v, see |^]). We will refer to the case v = 1 as 
nondegenerate, and to v > 1 as degenerate. 

5.2 The PRF reduction scheme 

We want now to consider the PRF corresponding to the linear part given by 
A. In the spirit of PRF, we should act on the NF (14) with Lie-Poincare 
transformations generated by homogeneous functions h m G Kei(Co)nV m . These 
will correspond to the action of vector fields of the form H m = aX m + (3Y m . 

We have then to consider C\\ this depends on the coefficients of the qua- 
dratic part Wi of the vector field W, which we write as W\ = aiXi + b{Yi. 

Under our assumption that v < fi, necessarily a\ — 0. In the nondegenerate 
case, bi / 0, while in the degenerate one, bi = 0. 

In the nondegenerate case we have W\ = &1Y1. We notice that [Yi,Xfc] = 
— Ffc + i and [Yi, Y k ] = 0; therefore Ker(A'Ji) = y. On the other hand, Ran(A / (i) 
also is given by y, and Ker(X^) = X. In this case we also have to consider 
higher order parts of W; the first step of the PRF procedure can eliminate all 
terms in Ran(A^i) and thus we will only consider terms in Ker(A^ ] f ). 

Let fi be as above, and let = a^X^ (all the Y k parts with k > 2 can be 
eliminated, as just recalled). Now M.^ is the restriction of £ M to Ker(.Mi) = 
Ker(£o) H Ker(£i): indeed the C m with 1 < m < /J, are zero and put no 
restriction. We have [X^^Yk] = fc^fc+p and thus Ker(A^ M ) = {0}: no further 
normalization is possible within the PRF scheme. 



X k := x 



1 3 x SW fc and Y k :=x k yd y eW k 



OO 




(14) 
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Thus the PRF is given in this case by 

oo 

W = Y + b 1 Y 1 +J2 "kXk , (15) 

fc=2 

where the hats on constants a k indicate that coefficients are in general not the 
same as those of the initial NF (14). 

We anticipate that the LRF reduction scheme can give a finite dimensional 
normal form for this case, see below. 

The previous discussion can easily be generalized to the degenerate case 
where 1 < v < /j, with /x and v defined above; notice that at least one of these 
has to exist and be finite, or the system would already be linear and thus trivial. 
Here the NF is 

\x— 1 oo 
k— v k—u 

Now C v (H k ) = b v ^Y v ,a k X k + (3 k Y k ] = -vb v a k Y u+k , and therefore we can 
eliminate all the Y u+k terms simply by choosing, with the same notation as 
before, a k = —b v+k /{vb v ); we cannot eliminate any of the X k terms. Thus, the 
PRF in the degenerate case is 

oo 

W = Y + b v Y v + "kX k . (17) 

k— p. 

Similarly to what happens for the nondegenerate case, the LRF procedure gives 
better results in this case. 

5.3 The LRF reduction scheme 

In the previous computations, we have followed the general PRF scheme for 
further normalizing the standard NF (14); this gave an infinite PRF in both the 
degenerate and nondegenerate cases. 

However one can take advantage of the specific Lie algebraic structure of 
Q = X®y, embodied in (13), to obtain a more drastical reduction: indeed, one 
can obtain a reduction to a finite normal form (the Lie renormalized form), as 
we now discuss. We use the same notation as in discussing the degenerate case 
above. 

We first operate a sequence of normalizations with generators h, a ^ = a k X k , 
which we choose so as to eliminate higher order X k terms, i.e. X k for k > fi 
(as we know, this is not possible for k = 2/i). Notice this will change not only 
the (coefficients of the) X k terms, but the (coefficients of the) Y k terms as well; 
however, no terms of degree k < v will be produced. 
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In this way, we arrive at a partially reduced form (the tilde indicates that 
the coefficients are not the same as the initial ones, but not yet final) 

oo 

W = Y + a^Xf, + a 2 ^X 2fi + J2~b k Y k . 

k—v 

Once this has been done, we pass to consider a second sequence of normal- 
izations with generators = [3 k Y k . As y is an ideal in Q, the X k terms are 
unaffected. On the other side, y is abelian, and so only the X M and X 2fl are 
actually active in these transformations: that is, we can only eliminate terms 
Y^+i and higher (it is clear by the commutation relations that these can always 
be eliminated). 

In this way we arrive at the LRF: this is a NF depending on [\i — v + 3) 
constants (recall we assumed [i > v\ of the form 

n 

W = Y + a^X^+^X^ + Y,bkY k . (18) 

k — V 

It is also clear by this discussion that actually b k — b k , a 2fi = a 2tl . Note that 
the number of constants (p — v + 3) agrees with that computed by Bruno JlT| . 

It should be stressed that this LRF is not a PRF, as can be checked by 
comparing this with (15) and (17) above, or comparing the definition of PRF. 

Indeed, with the notation employed in |2lJ, the spaces := (~)W k 

(k) 

with v < k < [i reduce to multiples of X k . Here we have therefore W k ; 
for v < k < /i, and thus the LRF cannot be a PRF. 

5.4 Explicit reduction (nondegenerate case) 

The reductions described in previous subsections can be explicitely performed; 
detailed computations are reported in [2^| , while here we just give results. We 
write the normal form W in the form (14) and consider further normalization 
up to order six. 

Let us first consider the PRF reduction. In the nondegenerate case we 
will take, for the sake of simplicity, [3 k — 0; we choose ol\ = —b 2 /bi, a 2 — 
{b\ - 6i& 3 )/6f, a 3 = -(261 - 3M 2 6 3 + &?&*)/&?, «4 = (96| - I8M2&3 + 36?6§ + 
86^6264 — 2b\b§) / (2b\) . In this way we obtain 

W<® = Y Q + 61 Yi + a 2 X 2 + [03 - a 2 b 2 /b x ] X 3 + 
+ [a 4 - 2a 3 6 2 /fei + a 2 b 2 2 /b\] X 4 + 
+ [a 5 - 3a 4 6 2 /fei + 4a 3 6^/&f + 

-a 3 b 3 /b 1 - 2a 2 b 2 b 3 /b\ + a 2 b 4 jb Y \ X 5 + 0(6) . 

Let us now pass to consider the LRF reduction, and perform detailed com- 
putations according to the LRF scheme for the nondegenerate case; that is, we 
first take care of the X k terms (up to k — 5), and then of the Y k ones. 
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With a transformation hi = a\X\, the W 3 term reads 

W 3 = [o 3 + a 2 ai] X 3 + [63 + 26 2 ai + haj] Y 3 . 

We disregard the Y 3 term and choose a\ so to eliminate the X 3 term, i.e. 
ai = — a 3 /a 2 . After computing the effect of this on higher order terms, we 
could perform a transformation with generator h 2 — a 2 X 2 . However, we know 
that there will be no way to eliminate the X4 term, so we set a 2 = 0. We 
perform a transformation with generator h 3 = a 3 X 3 . With this, the W5 term 
reads 

W 5 = [2a\/al - 3a 3 a 4 /a 2 + a 5 - a 2 a 3 ]X 5 + 

+ [a 3 b\/a 2 - Aa 3 b 2 /a 2 + %a\b 3 ja^ - Aa 3 b 4 /a 2 + 
+ b 5 - 2a 3 b 1 a 3 /a 2 + 2b 2 a 3 ]Y 5 . 

Again we only aim at eliminating the X5 term, and thus we choose a 3 = \{2a 3 — 
3a 2 a 3 a 4 + a 2 a5)/{a 2 )\- We will be satisfied with this order of normalization for 
the Xk terms, and take now care of the Yfc ones. 

We first operate a transformation with generator hi — PiY\; we get 

W 3 = [ajh/al - 2a 3 b 2 /a 2 + b 3 - a 2 /3i] Y 3 . 

By choosing j3\ = \{a 3 b\ — 2a 2 a 3 b 2 + a 2 b 3 )/(a 2 )] we eliminate this. We com- 
pute the effect on higher order term, and then consider a transformation with 
generator h 2 = [3 2 Y 2 ; with these, we have 

W 4 = [a 4 - a 3 /a 2 ] X 4 + 

+ (l/ol) [a 3 bi + 3a 2 a§6 2 - ia 2 a 3 (a 4 b 4 + a 2 b 3 ) + 
+4{a 5 h + a 2 (b 4 - 2a 2 (3 2 )) } Y 4 . 

We want to eliminate the Y 4 term, and thus we choose [3 2 = (l/2a 2 ) (a 3 b\ — 
3a 2 a 3 a 4 bi + a 2 a^b\ + 3a 2 a^b 2 — 3a 2 a3&3 + a 2 b 4 ). Again we take into account the 
effect of this on higher order terms, and pass to consider a transformation with 
generator h 4 — [3 4 Y 4 ; we get 

W 5 = [—2a 3 bi/a 2 + ha 3 a 4 b\/a\ - 2a 3 a 5 6i/a 2 - 2a|6 2 /a 2 i - 4a 3 a 4 b 2 /a 2 + 
+2a 5 b 2 /a 2 + 7a 3 b 3 /a\ - a 4 b 3 /a 2 - 4a 3 b 4 /a 2 + b 5 - 3a 2 /3 3 ] Y 5 

which can be eliminated by choosing (3 3 — — (l/(3a 2 )) (2a|fei — 5a 2 a^a 4 b 4 + 
2a 2 a 3 a 5 b 4 + 2a 2 a 3 b 2 + Aa 2 a 3 a 4 b 2 — 2a 2 a^b 2 — 7a 2 a 3 b 3 + a 2 a 4 b 3 + Aa 2 a 3 b 4 — a 2 b§). 

Summarizing, and having taken into account all higher order effects (up to 
order six), we have reached the LRF 

W = Y Q + hYt + a 2 X 2 + [b 2 - (a 3 h/a 2 )}Y 2 + [a 4 - {a 2 3 /a 2 )]X 4 + 0(6) 
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Appendix A. 

Bruno alternative definition and example. 

In his recent book Jl3| (and before this is HU), A.D. Bruno has dealt with 
PRFs. Unfortunately he reports a definition for PRFs which is not equivalent 
to the one contained in [^o[ ^lj (and reported above), so that his comments on 
PRFs can improve confusion rather than clarifying the issue, especially since 
the difference in definitions is not made clear (nor mentioned). 

The definition of PRFs given in Bruno's works (which I take from section 
V.22 of |l^| ) is as follows, once translated to the notation used in the present 
paper (the definition given in the first review jl^] is slightly different, but the 
difference is unessential here). 

We consider the vector power series ft] let F m = J2T=o f k > * ne operators V?*. 
are then defined as s S m (H) = {Wdj)F m — (F-> n dj)H; in the present notations, 
we have 

tf m = C k . (A.l) 

k=0 

Bruno considers then the adjoint of the operators $,„, and declares that 
PRFs have the property that 

n-i(Fk) = ; (A2) 

this is not equivalent to my original definition. In facts, this condition implies 
that each term Gk is in Ker(£+) for all the p < k; in my definition |20 , 21 this 



applies to Aip rather than £+, which obviously makes a substantial difference. 
Bruno considers then an example, i.e. the system 

x = x 3 , y = y(l + x + x 2 ) ; (A. 3) 

according to [jl3| the PRF for this is given by (0,2 7^ 0) 

x = a 2 x 3 +ax 5 , y = y(l + fix) (AA) 

with no higher order terms. In |L2| the PRF is instead claimed to be 

x = x 3 , y = y(l+x) , (A.5) 

again with no higher order term (the difference between these two expressions 
is never mentioned, nor explained, in both of jl2[ ft3|). In both cases, Bruno 
does not explain how these expressions are obtained. However, both (A. 4) and 
(A. 5) are obviously different from the real PRF, see the example discussed in 
detail in section 5. Bruno proves then that (A. 3) cannot be conjugated to (A. 4) 
or (A. 5); we have seen that these are not PRFs and thus his arguments do not 
concern PRFs (in the sense considered in [pol ^ll and here). 
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Appendix B. 

General Lie algebraic structure of normal forms 

A general approach employing the Lie algebraic structure of nonlinear vector 
fields resonant with a given linear one was proposed by Broer and Takens ^ , 
and further developed by other authors. The Broer procedure is always appli- 
cable, and is able produce a unique reduced normal form |^|, |^] . Unfortunately, 
implementing it in practice seems to be quite hard, and this approach has been 
applied only to system of very low dimension. On the other side, the central 
ideas put forward by Broer are at the basis of many approaches to the problem 
of further reduction of normal forms (in particular the one discussed in this 
note) and I will describe them the way I understand them. 

Given a matrix A, we denote by Q the Lie algebra of vector fields in normal 
form with respect to A. One can consider the descending central series (DCS) of 
Q [^7], |3^]. We recall this is defined by Go = G and Gk+i = [G, Gk]\ as well known 
the factor algebras Tk = Gk/Gk+i are abelian. It is then possible to eliminate 
terms by inner automorphisms of G (that is, by acting on G with vector fields in 
G) proceeding along Gk, i-e. filtering the Lie algebra G- This is relevant to our 
problem since Lie-Poincare transformations reduce to inner operations on G ■ 

This approach is, of course, particularly convenient when G is nilpotent; 
actually if we consider nonlinear vector fields resonant with a given A, i.e. G* '■= 
G\{G n V ), and Gk its DCS, we have Q% n V m = for m < k, and thus G* is 
nilpotent of infinite order. We can reduce to the more familiar case of algebras 
nilpotent of finite order TV by considering only the projections of Gk and Tk to 
the linear space C V spanned by nonlinear vector fields homogeneous of 
degree not higher than N. 

Let us now briefly describe the implementation of the BRF procedure. We 
consider the sequence of abelian factor algebras Tk = Gk+i/Gk, and denote by 
Xk the projection from G to IV 

We will produce a sequence of vector fields Wi , each being the result of 
the first j steps of further normalization. Let rjj be the operator on G defined by 

r)j(H) = e H W**^ e~ H ; let ttj be the operator of projection from G to the range 
of r]j. We can reduce the normal form W = wi " 1 as follows. 

As the first step, consider a £ Q and require that W = := r] (H) 

is such that 7r [xo 7 ?o(^^ ^)] = 0- This determines (non uniquely) and 
produces a wi X \ 

At further steps, we have the same setting; the "homological equations" on 
Lie algebras to be solved at each step will be 




(B.l) 



this determines H^h Each W* is then determined as 



Wi j+l) = Vi(H U) ) := cxp[H U) ] exp[Wi j) ] exp[-fT«] . 



(B.2) 
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